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Abstract 

We prove that the symbolic dynamical system generated by a purely substitu- 
tive Arnoux-Rauzy sequence is measurably conjugate to a toral translation. The 
proof is based on an explicit construction of a fundamental domain with a fractal 
boundary (a Rauzy fractal) for this toral translation. Such a fractal is obtained 
as the Hausdorff limit of a sequence of compact planar sets, where each set is the 
projection of a finite union of faces of unit cubes generated by a multidimensional 
version of the underlying Arnoux-Rauzy substitution. 

This paper is dedicated to the memory of G. Rauzy. 

1 Introduction 

Substitutions and self-similarity Self-similarity is a notion that has proved its rel- 
evance in the study of dynamical systems. Among the numerous ways of producing 
self-similar dynamical systems, substitutions play a prominent role. The term "substi- 
tution" encompasses symbolic replacement rules that act on words (generating symbolic 
dynamical systems), and more generally replacement rules acting on finite sets of compact 
subsets of R d (generating tiling spaces, the most celebrated example being the Penrose 
tiling). Substitutions prove to be very useful in many mathematical fields (combinatorics 
on words, number theory, harmonic analysis, ergodic theory), as well as in theoretical 
computer science and physics, such as illustrated by the bibliography of [PF02j . 

S-adic systems A natural generalization of substitutive symbolic dynamical systems 
arises when one considers the dynamical system obtained by iterating not one, but sev- 
eral substitutions taken from a finite set S of substitutions: such systems are said S-adic. 
For more details see, e.g., |Dur00| iLerlf j . These systems, which have hierarchical struc- 
ture, belong to the larger category of fusion tilings introduced in |PS11] , which also 
includes Bratteli-Vershik systems as well as multi-dimensional cut-and-stack transforma- 
tions. These hierarchies are not necessarily the same at each level (as it is the case for 
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purely substitutive systems), but there are only finitely many different possible structures 
at each level. 

Linearly recurrent systems form a widely studied class of systems that share many 
properties with substitutive dynamical systems. For instance, it has been proved that 
measure-theoretic and topological discrete spectrum are equivalent for linearly recurrent 
systems [CDHM03J, as it was already known before for dynamical systems arising from 
a single primitive substitution |Hos86] , see also |Sol07] . The notion of 5-adic system has 
been used notably for the description of linearly recurrent systems. For instance, Durand 
characterized linearly recurrent subshifts as the set of primitive proper S'-adic subshifts 
[DurOQllDur03j . 



Arnoux-Rauzy sequences Sturmian sequences are among the most classical examples 
of S-adic sequences and are a classical object of combinatorics on words and symbolic 
dynamics. They are the infinite sequences of two letters with factor complexity n + 1 
(i.e., they have exactly n + 1 factors of length n), and they correspond to natural codings 
of irrational rotations on the circle {MH40J. These sequences admit an S'-adic definition: 
they are the sequences that belong to the symbolic dynamical systems generated by 
iterating the two substitutions o\ : 1 H- 1, 2 H- 21, and a 2 : 1 H> 12, 2 H- 2. See [PF02I 
Chap. 6] and |Lot97| Chap. 2] for a detailed survey of their properties. 

Arnoux-Rauzy sequences were introduced in |AR91] in order to generalize Sturmian 
sequences to three-letter alphabets. They can be defined in an S'-adic way as the sequences 
that belong to the symbolic dynamical system generated by iterating the Arnoux-Rauzy 
substitutions o~i, a 2 , cr 3 given by 



a i ■ \ A , v AA if A -LA \ l = 




where each of the o~i occurs infinitely often in the iteration. These sequences have fac- 
tor complexity 2n + 1. It was conjectured that Arnoux-Rauzy sequences correspond to 
natural codings of translations on the two-dimensional torus (as in the Sturmian case 
with rotations on the circle), but this conjecture has been disproved in |CFZ00] . This 
family of sequences has indeed a complex and intriguing ergodic behaviour, as underlined 
in |CFM08j : some of them might even be measure-theoretically weak mixing. Arnoux- 
Rauzy sequences are also studied in word combinatorics, they belong indeed to the family 
of episturmian words (see the survey |GJ09j ). For more references about Arnoux-Rauzy 
sequences, see |BFZ05l ICC06] . 

In this article, we focus on purely substitutive Arnoux-Rauzy systems, that is, systems 
generated by an infinite purely periodic product of Arnoux-Rauzy substitutions (there are 
infinitely many such sequences). We focus on the spectrum of these dynamical systems, 
and our main theorem is the following (see Section [2] for precise definitions). 



Theorem 1.1. Let o be a finite product of the substitutions o\, a 2 , 03, where each substi- 
tution appears at least once. The symbolic dynamical system (X CT , S) generated by a has 
pure discrete spectrum. 



Our methods In order to prove Theorem 1.1 , we provide an explicit measure-theoretical 



isomorphism between the symbolic dynamical system generated by an Arnoux-Rauzy sub- 
stitution a and a two-dimensional translation of T 2 , with the fundamental domain of the 
translation being the Rauzy fractal of the substitution a. Rauzy fractals are compact 
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planar sets with fractal boundary that can be associated with any unimodular Pisot irre- 
ducible substitution. They first appeared in the work of Rauzy |Rau82j . who generalized 
the theory of interval exchange transformations by defining a domain exchange transfor- 
mation of three pieces in K 2 . These fractals prove to be very useful to investigate the 
properties of dynamical systems arising from Pisot substitutions (see Section 2.5 for more 
details). 

In this article, we exploit a combinatorial definition of Rauzy fractals introduced by 
Arnoux and Ito in [AI01J to obtain the discrete spectrum property. More precisely, we use 
methods inspired from [1094] : we prove that a particular sequence of approximations of 
the Rauzy fractal contains arbitrarily large balls, which is equivalent to the pure discrete 
spectrum property (see Section 4.2). 

Theorem has been announced in [AIQlj . in which it is stated that "we can check 
that the sequence of approximations contains a fixed open set", but no indication about a 
proof of this fact is given. Note that the "sequence of approximations" refers in [AlOl] to 
renormalizations of the sequence of approximations of the Rauzy fractal that we consider 
in the present paper. This explains why the authors of [AI01] are looking for a fixed open 
set contained in the approximations, whereas we prove that our approximations contain 



arbitrarily lar ge balls. Also, Theorem [Llj has recently been obtained by Barge, Stimac 
and Williams BSW11 using different methods based on the balanced pair algorithm and 
overlap coincidences. 

The present approach works by describing combinatorial and geometric properties of 
the Arnoux-Rauzy systems through the associated approximation sets, by introducing 
a finite graph (depicted in Figure [I] below) that governs some topological properties of 
these approximation sets. More generally, our approach is well-suited for handling infinite 
families of substitutions (as in Theorem 1.1) instead of single given substitution. Indeed, 
determining if a Pisot irreducible substitution has pure discrete spectrum is decidable 
|Sie04] [ST09, SS02 | lA"Lll] . but such decidability results are much more difficult to obtain 
for families of subs titutions , although some results have been obtained in this direction 
in [AkiOOl IBBK061 iBSWllj . Our approach can also be applied to S-adic systems (see 
Theorem 4.7). We thus plan to develop the methods described here to further 5-adic 
families of substitutions, such as Jacobi- Perron substitutions in [BB JSllj . 

Furthermore, the present approach provides geometric information on the associated 
Rauzy fractals (see Theorem 4.6 below and [BJS11] ). This has interesting consequences 
for exhibiting Markov partitions for toral automorphisms. Indeed, as a consequence of 
Theorem 1.1 , Rauzy fractals are known to provide explicit Markov partitions for the auto- 
morphism of the three-dimensional torus defined by the abelianization of the substitution 
a, for every finite product a of Arnoux-Rauzy substitutions (where each substitution a iy 
i = 1,2,3 appears at least once) [10931 IPra99| HR06J. Markov partitions are known to 
exist [Bow75j for such automorphisms, but no general explicit construction is known in 
dimension > 3. Note that the elements of the Markov partitions we provide associated 
with Arnoux-Rauzy substitutions are connected, by connectedness of the Rauzy fractals 
(see jBJSllj ). 



Outline of the paper Section [2] introduces basic notions and definitions. Specific 
combinatorial tools that allow the study of the topology of the patterns generated by 
Arnoux-Rauzy dual substitutions are described in Section [3j and in Section |4| we exploit 
the tools developed in the previous sections to prove our main results concerning the 
associated Rauzy fractals. 
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2 Definitions and preliminary results 



In this section, we introduce the basic objects of this article. We start with the classi- 
cal notion of substitutions on words, and then we define unit faces and discrete planes. 
This allows us to define dual substitutions, which act on discrete planes instead of uni- 
dimensional words. We then introduce Arnoux-Rauzy substitutions and Rauzy fractals, 
together with basic considerations on the spectrum of symbolic dynamical systems. 

2.1 Substitutions 

Let A = {1, . . . , d} be a finite alphabet, and A* the free monoid generated by A, that is, 
the set of finite words over A. 

Definition 2.1 (Substitution). A substitution is a non-erasing morphism of the free 
monoid A*, i.e., a function a : A* — > A* such that a{uv) = a{u)a{v) for all words 
u,v £ A*, and such that the image of each letter of A is non-empty. 

The incidence matrix M CT of a is the square matrix of size dxd defined by M CT = (m^), 
where ra^ is the number of occurrences of the letter i in cr(j). A substitution a is 
unimodular if det(M CT ) = ±1. The map t : A* — > N d ,w >-)■ (|w|i, \uj\2, ■ ■ ■ , \w\d) is called 
the abelianization map, where \w\i is the number of occurrences of the letter i in w. The 
map I is also referred to as the Parikh mapping. 

A substitution o is primitive if there exists an integer n (independent of the letters) 
such that a n (a) contains at least one occurrence of the letter b for all a, b E A. We will 
work with primitive substitutions having furthermore additional arithmetic properties. 
We recall that a Pisot number is a real algebraic integer greater than 1 whose conjugates 
have absolute value less than 1. 

Definition 2.2 (Pisot irreducible substitution). A unimodular substitution a is Pisot 
irreducible if the maximal eigenvalue of M CT is a Pisot number, and if the characteristic 
polynomial of M CT is irreducible. 

Let a be a unimodular Pisot irreducible substitution. We denote by Ug a positive 
eigenvector of associated with (3, and by a positive eigenvector of t M cr associated 
with (5. The existence of such vectors comes from the Perron- Frobenius theorem applied to 
M CT , because the matrix of a Pisot irreducible substitution is always primitive, as proved 
in [CSOl] . 

Definition 2.3 (Contracting plane). Let a be a unimodular Pisot irreducible substitution. 
The contracting plane V c associated with a is the plane of normal vector v^. We denote 
by 7r c : K 3 — > V c the projection of K. 3 on V c along ug. 

Symbolic dynamics of substitutions Let S stand for the (two-sided) shift on A z , 
i.e., S((u n ) ne z) = (u n +i)nez- We endow A z with the product topology of the discrete 
topology on each copy of A. Let us now see how to associate with a primitive substitution 
a symbolic dynamical system, defined as a closed shift invariant subset of A z . Let a be 
a primitive substitution over A, and let u G A z be such that a k (u) = u for some k > 1 
(such an infinite word exists by primitivity of a). Let 0(u) be the positive orbit closure 
of u under the action of S, i.e., the closure in A z of the set 0(u) = {S n (u) : n > 0}. By 
primitivity of a, 0{u) does not depend on u, so we can define 0{u) to be the symbolic 
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dynamical system generated by a, that we denote by X a . The system (X a , S) is minimal 
(every non-empty closed shift-invariant subset equals the whole set) and uniquely ergodic 
(there exists a unique shift-invariant probability measure on X a ). It consists of all the 
bi-infmite words that have the same language as u. See |QuelO| for more details. 

2.2 Discrete planes 

Let (e!,e2,e 3 ) denote the canonical basis of K 3 , and let R 3 denote the set of vectors of 
K 3 with positive entries. The following definition of a discrete plane was introduced in 
|IQ93t IIQ94] under the name of "stepped surface". The set of its vertices corresponds 



to the classic notion of "standard arithmetic discrete plane" in discrete geometry, as 
introduced in |Rev914 IAnd 03j . 

Definition 2.4 (Discrete plane). Let v G R 3 and S be the half-space of points x G R 3 
such that (x, v) < 0. The discrete plane T v of normal vector v is the boundary of the 
union of the closed unit cubes with integer coordinates that intersect S. 

A discrete plane can be seen as a union of unit faces [x, i]*, defined by 

[x,l]* = {x + Ae 2 + fie 3 : A, G [0, 1]} = ^ 
[x,2]* = {x + Xei + fj,e 3 : A, // G [0, 1]} = (X 
[x,3]* = {x + Aei +fie 2 : A,// G [0,1]} = <£> 

where j G {1,2,3} is the type and x G 1? is the position of [x, i]*, as illustrated below, 
ei ir^Ni e 2 



[(-1,1,0),!]* [(-3, 0,-1), 3]* 



We denote by x + / the translation of a face / by x, and we denote by x + P the set 
{x + / : / G P}, where P is a union of faces and x G Z 3 . 

Remark 2.5. Every discrete plane contains the "lower unit cube" defined by U = [0, 1]*U 
[0,2]* U [0,3]* =0. 

The following proposition is a useful characterization of discrete planes, that will often 
be used to prove that a given finite union of faces is not included in a given discrete plane. 

Proposition 2.6 f |ABI021 |A"BS04] ). The discrete plane T v is the union of faces [x,i]* 
satisfying < (x, v) < (e^, v). 

2.3 Dual substitutions 

We now introduce dual substitutions, which act on unit faces of K 3 instead of letters. 
This notion, which can be seen as a multidimensional counterpart of substitutions, was 
introduced by Ito and Ohtsuki |IQ93| IIQ94j and refined later by Arnoux and Ito [AlOlJ. 
A more general definition is given in iS.AlOlj. 

Definition 2.7 (Dual substitution). Let o : {1,2,3}* — > {1,2,3}* be a unimodular 
substitution. The dual substitution associated with a, denoted by E*(a), is defined by 

E*(a)([x,*]*) = M; 1 x+ U U [M; 1 ^),^, 

j'=1,2,3 s\cr(j)=pis 
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where £(s) G Z 3 is the abelianized vector of s. We extend this definition to any union of 
unit faces: E* 1 (a)(P 1 U P 2 ) = E*(cr)(Pi) U E*(a)(P 2 ). 



For examples of application of this definition, see Section |2.4[ Let us remark that 
E^(cr) is completely described by M CT and the images of the faces [0, 1]*, [0, 2]* and [0, 3]*, 
by linearity of the multiplication by M" 1 . Also, we have E*(cr o a') = E^cr') o E^(cr) for 
every unimodular substitutions a and a' (see |AI01j for more details). 

The proposition below describes the interplay between discrete planes and dual sub- 
stitutions. 

Proposition 2.8 ( |AI01t IFer06j ) . Let T v be a discrete plane and a be a unimodular 
substitution. The image ofT v by E^(cr) is the discrete plane Tt MaV . 

Furthermore, two distinct faces of T v have images by E^ (0) that do not have any unit 
face in common. 

Note that in the special case where T v = T Vfj is the discrete plane associated with the 
contracting plane of a (T Vfj is indeed a discrete plane, because has positive coordi- 
nates), Proposition 2.8 implies that Ej(cx)(r v/3 ) = Tt MijVi3 = T Vfj , which means that T V0 



is invariant under E^(cr). 



2.4 Arnoux-Rauzy substitutions 

Let 0i, 02, 03 be the Arnoux-Rauzy substitutions defined by 

f 1 H- 1 ( 1 i — y 12 [ 1 ^ 13 

0! : I 2 ^ 21 2 : < 2 ^ 2 3 : I 2 ^ 23 
[ 3 i— )■ 31 ( 3 ^ 32 [ 3 ^ 3 

and Si, S 2 , E 3 the dual substitutions defined by Ej = E*(0.j) for i G {1,2,3}. According 
to Definition |2.7[ we have 



Si([x. 


1]*) 


= M; 1 x + 


(0,0,0) 


,1]*U[(0,0,0),2]*U[(0,0,0),3]* 


Si([x. 


2]*) 


= M; 1 x + 


(1,0,0) 


,2]* 


Si([x. 


3]*) 


= M; 1 x + 


(1,0,0) 


,3]* 


S 2 ([x. 


1]*) 


= M; 1 x + 


(0,1,0) 


, 1]* 


S 2 ([x, 


2]*) 


= M; 1 x + 


(0,0,0) 


,1]*U[(0,0,0),2]*U[(0,0,0),3]* 


S 2 ([x. 


3]*) 


= M^x + 


(0,1,0) 


,3]* 


S 3 ([x. 


1]*) 


= M; 1 x + 


(0,0,1) 


, 1]* 


S 3 ([x, 


2]*) 


= M; 1 x + 


(0,0,1) 


,2]* 


S 3 ([x. 


3]*) 


= M~ 1 x + 


(0,0,0) 


,1]*U[(0,0,0),2]*U[(0,0,0),3]* 



which can be represented graphically as follows, where the black dot stands for x on the 
left hand side and for M o : 1 x on the right hand side. 

{ <£> ^ <> { <£> H> *>> l^>^® 

Proposition 2.9 QAI01J). Let a be a product of o\, 02 and 03 in which each substitution 
appears at least once. Then, is a Pisot irreducible substitution. 
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We will also need to "desubstitute", i.e., to compute preimages of faces. The following 
lemma gives the preimages of unit faces for the Arnoux-Rauzy dual substitutions. The 
proof, which is an easy enumeration of cases, can be found in Appendix [Aj By abuse of 
notation, E^" 1 ([x,i]*) stands for E A J" 1 ({[x, i}*}), i.e., for the union of faces [y,j]* such that 
EfcQy, j]*) contains the face [x, This notation will also be extended similarly to finite 
unions of faces in all that follows. 



Lemma 2.10 (Preimages of faces). Let x 

sr'([x,i]*) = 



(1) 



G 1? . We have 



x+y+z 

y 

z 

x+y+z 

y 



So 1 



X 

x+y+z 

z 

X 

x+y+z 

z 

X 

y 

x+yj+z 

y 

x+y+z 



u 



u 



u 



^x+y+z-1^ 



X 

y 

x+y+z— 1 



ifi 
ifi 

ifi 

Hi 

ifi 
ifi 



1 

2,3 



1,3 



2,3 



2.5 Rauzy fractals and spectrum 

Geometrical representations and spectra A fundamental question concerning sub- 
stitutive dynamical systems (X a , S) deals with the possibility of giving them a geometric 
representation. By a geometric representation, we mean a dynamical system of geometric 
nature that is topologically conjugate or measure-theoretically isomorphic to (X a , S). In 
particular, we are interested in conditions under which it is possible to give a geometric 
representation of (X a , S) as a translation on the torus (or more generally on a locally 
compact abelian group). 

This question can be reformulated in spectral terms. Indeed, one associates with the 
(uniquely ergodic) symbolic dynamical system (X a , S) (considered as a measure-theoretic 
dynamical system) the operator U acting on the Hilbert space C 2 (X, /i) (where \i stands 
for the invariant measure) defined as follows: 



U 



£ 2 (X,/i) 



£ 2 (X,/i) 
foT 



The eigenvalues and eigenf unctions of (X, T, /z) are defined as being those of the unitary 
operator U. The symbolic dynamical system (X a , S) is said to have pure discrete spec- 
trum if C 2 (X,fi) admits an Hilbert basis of eigenf unctions. According to Von Neumann's 
theorem, (X a , S) has pure discrete spectrum if and only if it is conjugate to a translation 
on a compact group. For more details, see |Wal82[ QuelOj . 

The question of the existence of a geometric representation can thus be rephrased as 
follows: which are the substitutions whose associated dynamical system has discrete spec- 
trum? Recall that measure-theoretic discrete spectrum and topological discrete spectrum 
are proved to be equivalent for symbolic dynamical systems associated with a primitive 
substitution jHos86j. 
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Rauzy fractals It is widely believed that every Pisot irreducible substitution has pure 
discrete spectrum; this is known as the Pisot conjecture (see |PF02l Chap. 7] and 
[BK06] for more details). A strategy for providing geometric representations (and tack- 
ling the Pisot conjecture) has been initiated by Rauzy in |Rau82j for the particular Tri- 
bonacci substitution a : 1 i-> 12, 2 i— > 13, 3 h> 1, where it is proved that (X a ,S) 
is measure-theoretically isomorphic to the two dimensional toral translation given by 
x x + (1/(3, 1//3 2 ) ((3 is the Pisot eigenvalue of M CT ). 

Rauzy fractals were first introduced in |Rau82j in the case of the Tribonacci substi- 
tution, and then in [Thu 89j in the case of the /3-numeration associated with the smallest 
Pisot number. They are good candidates to be fundamental domains for the toral transla- 
tion that is conjectured to be measure-theoretically isomorphic to (X a , S). Such fractals 
can more generally be associated with Pisot substitutions (see [BSQH IST09| IPF02] ). as 
well as with Pisot /3-shifts under the name of central tiles (see [Aki98|. IAki02j ). 

An effective combinatorial construction of Rauzy fractals, different from the original 
one, has been proposed in |AI01j . Let a be a unimodular Pisot irreducible substitution. 



Remark |Z5| and Proposition |Z8| enable us to iterate E*(a) onW= [0, 1]* U [0, 2]* U [0, 3]* 
in order to obtain an infinite increasing sequence of patterns (with respect to inclusion). 
It is then possible to project and renormalize the patterns ~E*(a) n (L() by applying M CT o7r c 



(see Definition 2.3), in order to obtain a sequence of subsets of the contracting plane V c 
that converges to a compact subset of V c . More precisely, for n > 0, let (V n ) n>0 be the 
sequence of approximations of the Rauzy fractal associated with a, defined for all n by 

V n = M a n o7T c oEl(a) n (U). 

Arnoux and Ito proved in |AI01] that (T> n ) n>0 is a convergent sequence in the metric space 
of compact subsets of V c together with the Hausdorff distance. This yields the following 
definition for the Rauzy fractal associated with a substitution. 

Definition 2.11 (Rauzy fractal). Let a be a unimodular Pisot irreducible substitution. 
The Rauzy fractal associated with a is the Hausdorff limit of the sequence of approxima- 
tions (D n ) n>0 . 

This definition allows us to give an equivalent formulation of the pure discrete spectrum 
property in terms of the approximations T> n of the Rauzy fractal. The "if" direction comes 
from |AI011HR05] and the "only if" direction has been proved in |BK06] (see also |BST10j ). 

Theorem 2.12. Let a be a unimodular Pisot irreducible substitution. The substitutive 
symbolic dynamical system (X a , S) has pure discrete spectrum if and only if the approxi- 
mation set T> n contains arbitrarily large balls, as n goes to infinity. 

Our aim is now to prove that the approximations T> n contain arbitrarily large balls 
when n tends to infinity, when a is a product of Arnoux-Rauzy substitutions. This will 
be the object of Sections [3] and 4.1 



3 Covering properties and annuli 



We now introduce the combinatorial tools that will be used in Section 4.1, in order to 
prove that that some increasing sequences of patterns contain arbitrarily large balls. 

A finite union of unit faces is called a pattern. If £ is a set of patterns, we denote 
by C = {v + P : P G C, v G Z 3 } the set of all the translations of patterns of C In the 
following, we will not make the distinction between a pattern and the set of its faces. 
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3.1 Coverings 



The covering property below (introduced in |I093j ) can be seen as an arc-connectedness 
type property with respect to a set of patterns C: every couple of faces of an £-covered 
pattern must be connected by a path of overlapping patterns of C. If C is well chosen, 
^-covering guarantees good topological properties for the sets V n , such as connectedness 
(see, e.g., |BJSllj for more details). 

Definition 3.1 (£-cover). Let £ be a set of patterns. A pattern P is C-covered if for all 
faces e, / G P, there exists (Qi, . . . , Q n ) G C such that: 

1. e G Qi and / G Q n ; 

2. Q k fl Qk+i contains at least one face, for all k G {1, . . . , n — 1}; 

3. QfcCP for all k G {1, . . . , n}. 

The following proposition gives a convenient combinatorial condition to prove that 
^-covering of the images of a pattern by a dual substitution. 

Proposition 3.2 ([I 093j ). Let P be an C-covered pattern, E be a dual substitution and C 
be a set of patterns such that E(Q) is C-covered for all Q G C. Then S(P) is C-covered. 



3.2 Strong coverings 



The notion of /^-covering defined in Section 3.1 is already enough to provide us with some 



of the topological properties which will be needed in Section 44 However, to ensure that 
these properties are preserved under the image of a dual substitution, we need a stronger 
property, namely strong /^-covering, that we define in this section and that is the main 
new ingredient of this article (this notion plays a key role in particular in the proof of 



(see also Remark 3.6) 



Lemma 4.3). Example 3.8 shows a concrete situation where strong £-covers are needed 



Before defining strong ^-covers, we introduce the set of patterns £ e d ge that consists of 
all the twelve edge-connected two-face patterns defined below (a pattern is edge- connected 
if it is £edge-covered) : 



£^ = [o,i]* u [(o,i,o),i]* 

= [0,1]* U [(0,0,1),!]* 



= [0,2]* U [(1,0,0), 2]* 

g = [0,2]* U [(0,0,1), 2]* 

= [0,3]* U [(1,0,0), 3]* 

<>> = [0,3]* U [(0,1,0), 3]* 



= [0,1]* U [0,2]* 

^ = [0,1]* U [(-1,1,0), 2]* 

= [0,1]* U [0,3]* 

g> = [0,1]* U [(-1,0,1), 3]* 

^> = [0,2]* U [0,3]* 

<fr = [0,2]* U [(0,-1,1), 3]*. 



Definition 3.3 (Strong £-cover). Let £ be a set of edge-connected patterns. A pattern 
P is strongly C-covered if 

1. P is ^-covered; 

2. for every pattern X G £ e d ge such that X C P, there exists a pattern Y G C such 
that X CY C P. 



The intuitive idea behind the notion of strong ^-covering is that every occurrence of 
a pattern of £ e d g e in P is required to be "completed within P" by a pattern of C. 
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Example 3.4. Define the patterns P = ^ and Q = < ^), and let £ = 
Both P and Q are ^-covered, but only Q is strongly £-covered because the two-face 
pattern X = Qy can be "completed" by ^> within Q but not within P. Another 



example of an £-covered pattern which is not strongly £-covered appears in Remark 3.11 



In the above definition of strong £-covers, we require X to be in £ e d ge because it is 
precisely what we need for the results of this article, due to technical reasons based on 
the particular form of the Arnoux-Rauzy substitutions, but the same kind of definition 
could be made with another set of patterns than £ edt 



Unfortunately, an analog of Proposition 3.2 does not hold in general for strong C- 



covers, but we will prove an ad hoc version of it in the particular case of Arnoux-Rauzy 



substitutions (see Proposition 3.10 below). 



3.3 Annuli 

The aim of this section is to define a notion of topological annulus in the setting of 
finite unions of unit faces. The main argument of this paper is based on this object: 



we will prove in Section 4.1 that, under some technical hypotheses, the image of an 



annulus by a dual substitution remains an annulus, thanks to an ad hoc generalization of 



Proposition 3.2 (namely, Proposition 3.10). This will enable us to deduce in Section 4.2 
that the approximations T> n of Rauzy fractals contain arbitrarily large balls when n tends 
to infinity. 

Definition 3.5 (Annulus). Let £ be a set of edge-connected patterns and T be a discrete 
plane. An C-annulus of a pattern P C V is a pattern ACT such that: 

1. P, A U P and T \ (A U P) are ^-covered; 

2. A is strongly ^-covered; 

3. A and P have no face in common; 



4. Pnr\ (pu A) = 0. 

Topological features of annuli are guarantied by Conditions ^ and Q. The intuitive 
meaning of Condition Q is that A completely "surrounds" P, in such that a way that 
P does not touch the "outside" region r \ (P U A). Conditions ([IJ and ^ provide a 
convenient combinatorial way to ensure that the topological features of Conditions ^ 
and Q are preserved under the image by Arnoux-Rauzy substitutions, as will be proved 
in Section |47Tl 



Remark 3.6. Contrarily to what is claimed in [1094] , the image by a dual substitu- 
tion of a pattern that satisfies Conditions ( TJ) , (|3]) and Q has no reason to satisfy again 

3. 8| below. Nevertheless, the notion of strong 



these conditions, as illustrated by Example 
covering with respect to a well-chosen set of patterns will alllow us in the case of the 
Arnoux-Rauzy substitutions to prove that the image of an annulus remains an annulus. 



Example 3.7. Let £ar be the set of patterns defined in Section [374] , and let 
p i u M = and = 





where Pi = P2 = U is shown in dark gray and Ax, A2 are in light gray. (The underlying 
discrete plane T is not shown.) The pattern A\ is not an £AR-annulus of Pi, because 
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Condition [4] fails (the lower vertex of U intersects V \ (Pi U Ai)). However, A 2 is an £ar- 
annulus of P^. Conditions ([!]) and ([2j) are easy to check and Conditions (|3| and Q clearly 
hold. 



Example 3.8. Let P U A U / be the pattern shown below. 



PUAUf 




£2 £9 




(1) 



(P is shown in dark gray, A in light gray and / in white.) The pattern A is not an 
^AR-annulus of P, because it is not strongly £AR-covered. The strong £AR-covering fails 
because of the two-face pattern of A shown below in white: 




We can observe that £2^2(^4) does not satisfy the topological requirements needed to be 
an £AR- ann ulus of S 2 S 2 (P). However, if we "complete" the faulty pattern <<> by < ^> 
in P U A using the single white face /, then A becomes an £AR-annulus of P, as well as 
its image by £2^2, which also becomes "completed" by the faces in white as shown in 
Equation ([!]) above. 

3.4 Covering properties for Arnoux-Rauzy substitutions 

Let £ar be the set of patterns consisting of (X] , £0 , <£] , ^f" , ^> , < ^ £ £ e d ge and of 

= [0, 1]* U [(0, 1, 0), 1]* U [0, 3]* ^ = [0, 1]* U [(0, 0, 1), 1]* U [0, 2]* 

^> = [0, 2]* U [(1, 0, 0), 2]* U [0, 3]* fa = [0, 2]* U [(0, 0, 1), 2]* U [0, 1]* 

= [0, 3]* U [(1, 0, 0), 3]* U [0, 2]* <^ = [0, 3]* U [(0, 1, 0), 3]* U [0, 1]*. 

The goal of the remaining of this section is to show that the set of patterns £ar 
behaves well with respect to Arnoux-Rauzy dual substitutions: Z^R-covering is preserved 



by applying Si, S 2 or £3 (Proposition 3.9), and strong £AR-covering is also preserved 



(Proposition 3.10), under a pattern-avoidance condition that is always verified in the 



cases that we will encounter in Section 4.1 (Lemma 3.12). 

Proposition 3.9. Let P be an C^R-covered pattern. Then the pattern Sj(P) is Car- 
covered for every i G {1, 2, 3}. 



The proof (which is detailed in |BJS11] ) relies on Proposition 3.2 Sj(P) is £ar- 
covered for all P G £ar, which can be checked by inspection. The following stronger 



property will be needed in Section |4.1[ to prove the annulus property for Arnoux-Rauzy 
substitutions. 

Proposition 3.10. Let P be a strongly C^R-covered pattern that does not contain any 
translate of one of the three patterns 

^ = [(0,0,0),1]*U[(0,1,1),1]*, 
0d = [(0,0,0),2]*U[(1,0,1),2]*, 
^ = [(0,0,0),3]*U[(1,1,0),3]*. 

Then £«(-P) is strongly C^R-covered for every i G {1, 2, 3}. 
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Proof. Let E = Ej for some i G {1,2,3}. By Proposition 3.9, E(P) is Z^R-covered. Now 



let X C E(P) be an element of C e d ge - If X is a translation of one of the six patterns 
^ ' <^ ' ^ ' ^> ' ^ ' then the strong covering condition is trivially verified because X is 
itself in £ar- 

It remains to treat the case of the six remaining possibilities for X up to translation. 
We have X C S(P), so it is sufficient to check that for every pattern Xq C P such that 
X C E(Xo), there exists K G £ar such that X C K C E(Jfo). Moreover, since X is a 
two-face pattern, we can restrict to the case where Xq consists of two faces only, which 
leaves a finite number of cases to check for X . 

Suppose that X = t^. The following table displays all the possibilities for Xq. It has 
been obtained by applying Lemma 2.10 to the two faces of X. The first column of this 
table provides the index i for which E = Ej. 



% 


X such that X C E«(X ) 


E,(X ) 


2 




ft] 


2 






3 


<C^= [0,3]* U [(0,1,0), 1]* 




3 







In the first and third rows, there is indeed a pattern Y G Car such that X C Y C E(Jfo), 
namely F = , which settles the case of these two rows. The case of the fourth row is 
settled directly because we have assumed that P does not contain the pattern hf^l. For 



the second row, since P is itself strongly Z^R-covered, X is contained in some Y G Car 
(here, <j>^)- Hence, S 2 (lo) C E 2 (P), with E 2 (lo) G Car, so the property holds. 

The five remaining cases for X can be dealt with in exactly the same way as above: 
in each of the five cases, there are four preimages to check, which are symmetrical copies 
of the ones in the table above. □ 



Remark 3.11. The condition of avoidance of the three patterns in Proposition 3.10 
cannot be dropped, as illustrated by the following example: 




3> 




Indeed, the left pattern P is strongly £AR-covered and contains but S 2 (P) (depicted 
on the right) is not strongly ^AR-covered because of the two-face edge-connected pattern 
shown in dark gray, which cannot be included in any pattern Y G Car such that Y C P. 

We can observe that the pattern P above has a "hole". If we "fill" this hole with 
the corresponding unit face, P would no longer be strongly ^AR-covered. It would be 
tempting to try to put some topological restrictions on P to get rid of the assumption of 



avoidance of the three patterns in Proposition 3.10 but we will actually need to apply 



Proposition 3.10 to some patterns that are homeomorphic to topological annuli (as P is), 
in the proof of Lemma 4.3| 



In order to be able to apply Proposition 3.10 to patterns obtained by iterating Ei, 
E 2 and E 3 from U, we need the following lemma, which states that none of the three 
above-mentioned forbidden patterns appears in the iterations. 

Lemma 3.12. Let E be an arbitrary product of Ei ; E 2 and E3. The discrete plane 
^(r (1,1,1)) contains no translate of any of the three patterns 00, <<^. 
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Proof. We prove the result by induction on the size of the product E. First, we observe 
that r^i i) does not contain any of the three patterns in question: consider for instance 
a translate of the pattern [^] = [x, 1]* U [x + (0, 1, 1), 1]*, for x G Z 3 . We have (x + 
(0,1,1), (1,1,1)) = (x, (1,1,1)) +2 > (ei, (1,1,1)), which contradicts Proposition |2"!6 
The same reasoning holds for the two remaining patterns 00 and g for the initialization 
of the induction. 

In the table below, we have listed all the preimages of the three patterns (in light 
gray) obtained by applying Lemma |2.10| together with their only possible "completion" 
within a discrete plane (in dark gray). This completion is deduced from the arithmetic 



description of discrete planes provided by Proposition 2.6 



i 




% 




% 




2 




1 




i 




2 




1 




i 




3 




3 




2 




3 




3 


% 


2 





Let us work by contradiction and assume that a translate of one of the three patterns 
appears in E^i^i)) but not in E'(r (1,1,1)), where E = Ej o E', for some i e {1,2,3}. 
The above table, together with the injectivity of Ej (see Proposition 2.8), implies that a 
translate of one of the three patterns must appear in E'(P( 1)1 >1 )), which yields the desired 
contradiction. □ 



4 Main results 

We now use the results of Section [3] to prove the annulus property for Arnoux-Rauzy 
substitutions, from which we will deduce our main theorem and some other dynamical 
consequences. 

4.1 The annulus property for Arnoux-Rauzy substitutions 

The aim of this section is to prove the following theorem. The proof (by induction) relies 



on Lemma 4.2 (the base case) and on Lemma 4.3 (the induction step). 

Theorem 4.1. Let E be a product o/S 1; E 2 and E 3 in which each substitution appears 
at least once. Then, E n+2 (W) \ E n (W) is an CAR-annulus °/E n (W) in the discrete plane 
E n+2 (r(i i i i i)), for all n G N. 



Proof. We prove the result by induction on n. The case n = is settled by Lemma 4.2 



Now, assume that the induction property holds for some n e N. According to Proposition 



2.8 



E n+2 (W) is contained in the discrete plane E n+2 (r(i i i i i)). Hence, Lemma 3.12 allows us 



to apply Lemma 4.3 to deduce that E n+3 (W) \ E" +1 (W) is an £ A R-annulus ofE^+^ZY). □ 



Lemma 4.2. Let E be a product o/Ei, E 2 and E3 in which each substitution appears at 
least once. Then, S 2 (W) \ U is an Car- annulus ofU. 
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Proof. Let w — i\ ■ ■ -i n G {1, 2, 3} n such that E = E^ • • ■ E in . This proof is a case study, 
which we formalize as the study of the directed graph in Figure [TJ The vertices of this 
finite graph are patterns. This graph is such that for every edge P Q we have 
Q C Ej(P), and every vertex has three outgoing edges with distinct labels 1, 2 and 3 (for 
the sake of clarity, loops are not drawn in Figure [I]). It can be checked that one of the 
six bottom vertices is always reached when following the path of labelled edges given by 
w 2 = ww starting at the top vertex IA. 

The result follows because each bottom vertex is a pattern that consists of U and an 
^AR-ammlus of U, and because each outgoing edge of a bottom vertex points to a bottom 
vertex (these edges are not drawn in Figure [T]) , so applying E twice guarantees us to reach 
an £AR-annulus of U. 




Figure 1: Each vertex has three outgoing edges with distinct labels 1, 2 or 3. The 
18 outgoing edges of the bottom vertices are not drawn, but they all point to 
bottom vertices. The 9 other edges which are not drawn are loops. 



Let us note that the existence of a finite graph such as the one in Figure [TJ that provides 
the description of the possible annuli surrounding the unit cube is the key ingredient of 
the proof. There is no reason a priori for such a finite graph to exist for a given set of 
substitutions. 



Lemma 4.3. Let Y be a discrete plane that does not contain any of the patterns 
Ofl, 1>>- Let A C T be an CxK-o,nnulus of a pattern P C Y, and let E = Ej for some 
i G {1,2,3}. Then E(A) is an Cp^-annulus o/E(P). 

Proof. We must prove the following: 

1. E(P), S(v4) U E(P) and Y \ (E(A) U E(P)) are £ A R-covered; 

2. E(A) is strongly £AR-covered; 

3. S(A) and E(P) have no face in common; 



4. E(P) n E(r) \ (E(P) U E(A)) = 0. 
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Conditions ([I]) and ^ hold by Propositions |3. 2| a nd 2.8 respectively, and ^ holds thanks 
to Proposition 13.10 It remains to prove that (|4|) also holds. 



We will work by contradiction and we suppose that Q does not hold. This implies 
that there exist faces 

feP, g er\(AuP), feS(/), 

such that /' and g' have a nonempty intersection. Note that / U g is disconnected because 



P and r \ (P U A) have empty intersection by hypothesis. (See Figure |2 

0° 






Figure 2: AU P (left) and T,(A) U E(P) (middle or right), under the assumption 
that Condition (lil) does not hold. 



The strategy of the proof is as follows: we check all the possible patterns fUg and /'Up' 
as above, and for each case we derive a contradiction. This can be done by inspection of 
a finite number of cases. Indeed, there are 36 possibilities for /' U g' up to translation 
(the number of connected two-face patterns that share a vertex or an edge), and each of 
these patterns has a finite number of two-face preimages given by Lemma |2. 10 

The first such case is the following: 




[0,1]*U[(1,-1,1),1]* 
[0, 2]* U [(1,-1,1),!]*. 



By Proposition |2.6[ the pattern /Up can only be "completed" in the following way 
within the discrete plane T: 

fUgUX = 

where X = [0, 2]* U [(0, 0, 1), 2]* is shown in dark gray. Moreover, we have ICi because 
Condition Q for A and P would fail otherwise. However, this situation is impossible: 
X G £ c d g c but there cannot exist a pattern Y e Car such that X C Y C A because Y 
would overlap with / or g, which both are not in A. This implies that A is not strongly 
^AR-connected, which is a contradiction. 

The same reasoning applies to all the other possible cases for /' and g', as shown by 
the table of Figure |3} (By assumption, we have not considered the (problematic) patterns 



4.2 Dynamical consequences 



Our main theorems now follow easily from Theorem 4.1 



Theorem 4.4. Let a be a finite product of the Arnoux-Rauzy substitutions o\, a%, 03 in 
which each substitution appears at least once, and let E = E*(er). Then, for all R > 0, 
there exists n > such that the ball of radius R centered at the origin is contained in 
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/'US' 


/Ug 
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Figure 3: Disconnected preimages of connected two-face patterns by Si (left), S 2 
(middle) and S 3 (right). The only possible completion of fi)g is in dark gray, and 
/ U g, f U g' are in light gray. 



Proof. Let n > and T v be the discrete plane associated with the contracting plane of 



cr (we have £ n (£Y) C T v for all n > 0). By Theorem 4.1, the pattern T? n (U) consists of 
U and n disjoint £AR-annuli Ai, . . . , A n such that A i is an £AR-annulus of U and Ak is 
an ^AR-annulus of U U A\ U • • • Ak-i for all 2 < k < n, in the discrete plane r v/r 

Let us define the combinatorial radius of S 2n (W) as the smallest edge-connected dis- 
crete path of unit faces from a face of U to a face of T v \ S 2n (W). Condition (Jll of 



Definition 3.5 ensures that T? n {U) is simply connected (because its complement IsTar- 
covered and the patterns of £ar are edge-connected), so Conditions (|3| and ^ ensure 
that each annulus adds at least 1 to the combinatorial radius of T? n {U) (again, because 
the patterns of £ar are edge-connected). 

It follows that the combinatorial radius of T? n {U) is at least n, which implies that 
7r c (£ n (Z/)) contains arbitrarily large balls when n goes to infinity. □ 



We can now apply Theorem 2.12 to obtain Theorem |1.1 



Corollary 4.5. Let a be a finite product of the Arnoux-Rauzy substitutions o~\, a 2 and 
03 in which each substitution appears at least once. Then, the symbolic dynamical system 
generated by a is measurably conjugate to a toral translation, i.e., it has pure discrete 
spectrum. 

We also derive further topological properties of the associated Rauzy fractal. 

Theorem 4.6. Let a be a finite product of the Arnoux-Rauzy substitutions 0\, 02 and 03 
in which each substitution appears at least once. The Rauzy fractal associated with o is 
connected and is an inner point. 

The connectedness has been established in |BJSllj . and the fact that is an inner 



point comes from Theorem 4.4 together with [ST 09J. 



More generally, we have the following "S-adic version" of Theorem 4.4 
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Theorem 4.7. Let u be an Arnoux-Rauzy sequence defined as 



u = lim a io o ■ •■ o a in , 

n— s-oo 



where the sequence (i n ) e {1, 2, 3} N ; the substitutions Oi (i = 1, 2, 3) are the Arnoux-Rauzy 
substitutions, and the sequence (i n ) ne N takes each value in {1,2,3} an infinite number of 
times. Then, for all R > 0, there exists n > such that the ball of radius R centered at 
the origin is contained in 

7T c (E*K)o...oE* 1 ( ( T ? ;J(W)). 



Proof. The proof follows the same lines as the proof of Theorem 4.4 Let (nk)ke_N be 
an increasing sequence of integers such that the word i nk +i • • • in k+1 contains each of the 
letters 1,2,3 for every k. The role played by the substitution a will be played here by 
°i„ fc +i • ■ ■ ° °in fe+1 » and the role played by S will be played by E* (cr infe+i ) o ■ ■ ■ o E* (a ink+1 ) . 
Let 

s* = E;((7 j0 )o...oE;(aO- 



One proves exactly as for Theorem 4.1| that Sfc +2 (W) \ Efe(ZY) is an £AR-annulus of Sfc(W) 
in the discrete plane S fc+2 (r(i,i,i)), for all k G N. Indeed, by definition of the sequence 



(rifc)fc gN , in the proof of Lemma 4.2, any path labelled by i nk +i • • • ^n fc+2 starting at the top 



vertex U reaches one of the six bottom vertices. □ 

Such a result can be seen as a generation method (in the vein of |Fer08t IBLPPli] ) 
for a discrete plane whose normal vector is in some sense dual to the vector of letter 
frequencies associated with the infinite word u. It can also be considered as a first step 
towards the definition of Rauzy fractals and Rauzy tilings in the S"-adic framework. 



A Proof of Lemma 2.10 



Proof. We treat the case of Si only, the other cases being analogous. The images of the 
faces by Si are given by 



i — y 



(T*) + 

/ x-y-z\ 
{ I ) + 



x-y-z 

y 



U 



U 



A face of type 1 appears only once in the image of a face of type 1. If [(x, y, z), 1]* C 
Ei([(x ,yo,zo),l]*), then 



x 
V 



Hence 



xo - yo 



x + y + z 

y 

z 



S^aOr,^),!]* 



x+y+z 

V 

z 



A face of type 2 appears once in the image of a face of type 1 and once in the image of a 
face of type 2. If [(x,y,z),2\* C Si([(x , y , zo), 1]*), then 



x 

y 



%0-yo- z 
yo 

zo 



x + y + z 

y 

z 
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If [{x,y,z),2}* C Sx([(a;o,?/o,2o),2]*), then 





Hence 



Er 1 ^,^),^ 



2/ 



u 



x+y+z—1 
V 



A face of type 3 appears once in the image of a face of type 1 and once in the image of a 
face of type 3. In the same way as for a face of type 2, we get 



E-\[(x,y,z),3}*) 



( x+y+z 

[ y 

\ z 



)>1 



U 



^ X--!:-- . - 



z / 
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